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Abstract 

We give arguments that in the 1+1 dimensional abelian Higgs 
model the classical approximation can be good for the leading high 
temperature behavior of real time processes. The Chern-Simons dif- 
fusion rate ('sphaleron rate') is studied numerically in this approxi- 
mation. New results at high temperature show a T 2 / 3 behavior of the 
rate at sufficiently small lattice spacing. 
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1 Introduction 



In theories of baryogenesis the rate of sphaleron transitions plays an impor- 
tant role [[l], § . This rate is hard to calculate analytically, especially at high 
temperatures. Numerical simulations in the quantum theory have to face the 
problem of a complex 'Boltzmann factor' associated with real time correla- 
tion functions at finite temperature. For this reason a classical approximation 
was proposed some time ago ||, which has been tested in the abelian Higgs 
model in 1+1 dimensions [4-9]. At low temperatures the rate was found 
@, [7], H to agree with a semiclassical analytical calculation |I| , while at high 
temperatures the results have been somewhat confusing: ref. || argued for 
a T 2 / 3 behavior, while we argued |J for T 2 behavior. We found in fact T 2 
behavior at nonzero lattice spacing which turned into T-independence in the 
limit of zero lattice spacing. In the mean time simulations of the physically 
relevant SU(2) models in 3+1 dimensions have been improving [11-16] and 
the understanding of the classical approximation has been steadily increasing 
[17-24]. It therefore appropriate to return to the abelian Higgs model and 
try to assess the situation in 1+1 dimensions. We give a more full account 
of our earlier work briefly reported in ||. 

In sect. 2 we introduce the model and the rate of fermion number vio- 
lation and we discuss some of their properties. The classical approximation 
is introduced in sect. 3 and studied in perturbation theory. Sect. 4 and 5 
give details on the numerical method, sects. 6 and 7 contain results on the 
rate in the regimes of low and high temperatures and sect. 8 contains our 
conclusions. Details on the simulation algorithm using 'Kramers equation' 
are in the appendix. 



2 Abelian-Higgs model 

The classical abelian-Higgs model is given by the action 

1 



S = - / d 2 x 
with F, 



4g 



2 F^F^ + (D^)*Dy + /i 2 |0| 2 + A|0| 4 



[iv - d^Au - d v Afj, and F)^ = (<9 M - iA^cf) (our metric is given by 
goo = —Qn = — 1)- We recall that A^, \f~\ and g have dimension of mass, 
while <p is dimensionless in 1+1 dimensions. The role of a dimensionless 
coupling is played by A/|/x 2 |, at given £ = g 2 /A. 
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The classical model has a Higgs phase for ji 2 < 0, where gets a ground 



state expectation value |0| 



v/y/2, v 2 = — /i 2 /A, with masses mh = 1\v^ 



2„,2 



9 v 



In the euclidean version of the quantum theory there is a 
Berezinsky-Kosterlitz-Thouless phase transition in the scalar model at g = 
which is turned into a crossover at nonzero gauge coupling |25], p6fl . For a re- 
cent lattice study, see [27]. At nonzero temperature there is only a crossover 
in any case because then the system is infinite in one dimension only 

A toy model for the electroweak theory is obtained by coupling to fermions, 



S f 



d 2 x 



^(d, + i-A^ 5 )iP + yi>{(j>P R + <j>*P L )ip 



(2) 



where Pl,r — (1 T 7s)/2 are the chiral projectors. In the quantum theory 
the fermion current is anomalous, 



d^vf-ty = -q , 
with q the topological charge density (eoi = +1" 



(3) 



—e^F, 

4:71 



where the Chern-Simons current 

C » = ^ V A V . 
As a consequence a change in Chern-Simons number 



C 



dx C° 



o 



l r L 
— / dx A 1 
2ir Jo 



is accompanied by a change in fermion number 

rL 







dx ipi^°ip 



such that 

F(t)-F(0) = -[C(t)-C(0)} 



/ dt' / dx q(x, t') 
Jo Jo 



(4) 



(5) 



(6) 



(7) 



(8) 



We take space to be a circle with circumference L. 
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The sphaleron rate V can be identified from the diffusion of Chern-Simons 
number [|28|j , 



r 



C(0)f 



t 



oo. 



(9) 



In the following we shall neglect the influence of the fermions on the Bose 
fields. Our task is to evaluate the above real time correlation function. 

Transitions with AC of order 1 have to cope with the sphaleron barrier of 
energy E s = (2/3)v 2 m ( f ) , for which C = 1/2 p9 |. The sphaleron configuration 
can be written in the form 

TT V ms 7T 7T 

A = 0, A i = j, 0= y|tanh(— ?/)exp(2-y-z-), 

where y = x — L/2 G (— L/2, L/2). At relatively low temperatures <C 
T ^ E s the rate has been calculated semiclassically with the result 
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miL 




/(0 



(2tt) 3 

2i, 



>+l) 



s — - 



T(a + s + l)T{a- 
T(a + l)r(a) 

1 - 



nl/2 



1 



(10) 
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At high temperatures the rate is difficult to calculate analytically. On 
dimensional grounds it seems reasonable to expect that for temperatures T 
larger than any mass scale, in particular T ^> E s , 



(13) 



with dimensionless k. This formula is analogous to dimensional T/L 3 = 
k(o.wTY in 3+1 dimensions, with v~ 2 playing the role of the electroweak 
coupling aw- However, the behavior (|13|) is not obvious, since the original 
couplings A and g 2 are dimensionful. At high temperature the dimensional 
reduction approximation leads to temperature appearing in the combinations 
AT and g 2 T. Since these have engineering dimension three in mass units, we 
may expect the behavior 



^ = £(0(AT) 2 / 3 , 



(14) 



instead of ([13]). 
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3 Classical approximation 



The rate T defined by @ is a nonperturbative quantity and one approach to 
its computation is by numerical simulation. This is well-known to be difficult 
since the effective Boltzmann factor in real time processes is complex: 

Tr e~ H l T \e iHt C(0) e~ lHt - C(0)] 2 
(iC(t) ~ C(0)] 2 ) = L Tre J /r ^ (15) 

To cope with the complex weights a classical approximation has been in- 
troduced 0, ||] in which the quantum mechanical expectation value ( p3| ) is 
replaced by a classical expression, 

([C(t) - C(0)] ) = JD<pDire-B<*Wr " (16) 

Here ip and 7r denote generic canonical variables and ip(t) and 7r(t) are solu- 
tions of the classical Hamilton equations with an effective hamiltonian H e g 
and initial conditions <f{0) = <f, vr(0) = n. This approximation is used at 
high temperature, where the important low momentum modes have high 
occupation numbers. 

The classical approximation has been studied in perturbation theory in 
3+1 dimensions [17-21]. For static quantities the situation is well under- 
stood, as this corresponds to dimensional reduction in the imaginary time 
formalism, of a 4D theory to a 3D theory which can be renormalized in the 
usual way. For time dependent quantities there are divergent effects in clas- 
sical gauge theories (related to the physics of Landau damping) which are 
non-local and cannot be absorbed by local counterterms. In the renormal- 
izable quantum theory such effects correspond to loop momenta of order of 
the temperature - the well known hard thermal loop effects. 

In 1+1 dimensions the UV-divergencies are less severe. We shall now 
study the classical approximation for the abelian Higgs model in perturba- 
tion theory, using the imaginary time formalism and making an analytic 
continuation to real time at a suitable point. For simplicity we choose the 
classically symmetric phase fi 2 > for the starting point of perturbation 
theory. As a consequence the bare A-mass is zero, which causes infrared 
divergences in the diagrams. Such infrared divergences are presumably an 
artefact of perturbation theory and we assume they are cured automatically 
in a correct nonperturbative treatment. Here we shall use brute force and 
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introduce an infrared regulator mass m A . This crude method breaks gauge 
invariance but it is adequate for our purpose of elucidating the classical ap- 
proximation. 




Figure 1: Higgs and gauge field selfenergy diagrams. 

selfenergy given by the 



Consider the gauge boson contribution to the 
appropriate diagrams in Fig. |]: 

{k + 2p) fl {k + 2p) u 



E(p) 



D-U 



\J_.V 



d u -'k 



(17) 



(18) 



{2Ti) D - 11 

where D (— > 2) is the number of spacetime dimensions, kp = n2irT are 
the Matsubara frequencies (n = 0, ±1, ±2, . . .), and the vector propagator is 
given by 



D^(k) 



8^ + {a 



k^k v 



m l A + k 2 



(19) 



For illustration we first choose the Feynman gauge a = 1 and follow a stan- 
dard method f30| , which consists of carrying out the sum over Matsubara 
modes in the loop diagrams and subsequently analytically continuing the 
time component pn of the external momentum to real frequencies p°. Writ- 
ing (k + 2p) 2 = 2[m 2 + (k+p) 2 ] — [m 2 A + k 2 ] — (2m 2 —m 2 A — 2p 2 ), to conveniently 
cancel the time component of the loop momentum in the numerator against 
the denominators, we get 




x 



2uj + 
(1 +n + ri)(u + u') 



l + 2n . l + 2n' . 2ml - m\ - 2p 2 + 2p 2 



2uj' + 
(n — n'){ui 



2uouo' 



UJ 



(to + OJ') 



A2 



Po 



u') 2 -pi 



(20) 
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where 



11 = T) f = 

e"/ T -V e"'/ T -V 



10 



m 2 A + k 2 , J = Jm 2 + (k + p) 2 . (21) 



The external real-time frequency p° = —p Q is supposed to have a small pos- 
itive imaginary part, p° — > p° + ie, which corresponds to retarded boundary 
conditions. 

The zero temperature contribution, obtained by letting n,n' —> 0, has 
a logarithmic UV-divergence which can be canceled by a mass counterterm. 
The temperature dependent contribution is finite as usual. We now add the 
mass counterterm to E, assume dimensional renormalization in the MS-bar 
scheme, and denote the resulting selfenergy by E. We are interested in its 
high temperature behavior when T ^> m^, m^, |p|, p°. 

For high temperatures it turns out that the leading contribution in a high 
temperature expansion is obtained by substituting the leading behavior of 
the Bose distribution functions n and n', 

1 T 1 , T 

2 + "-^ 2 + "^' (22) 
The resulting expression can be written as 

D-2 1 2m\ -m\- 2p 2 + 2p 2 




uo 2 u 12 2lu 2 u' 2 

(23) 



(to + to') 2 (uj-uj') 2 



(uj + uj') 2 -pl (uj-uj') 2 -pl 

This is finite for D = 2, which justifies the replacement (^) under the inte- 
gral. This contribution is actually identical to the classical selfenergy which 
follows from (|16D, hence the subscript c. Intuitively this is plausible, since 
large occupation numbers n, n' imply classical behavior, although it is of 
course not obvious that solving the classical equations of motion and aver- 
aging over initial conditions according to ( |TED leads to the same correlation 



functions as making the high temperature approximation (0). See |^0[ for 
more details in scalar field theory. 

For p° = the classical selfenergy E c goes over into the selfenergy of 
the naively dimensionally reduced theory, the one dimensional field theory 
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obtained by keeping only the zero mode (n = 0) in In this case S is 

i it 

the selfenergy of the zero mode J dxD<fi(x, xr>). 

The difference between the classical and quantum selfenergy is subleading 
for large T and it can be expressed as 



9_ 

2tt 



in|) -i- ocr 1 - 



(24) 



-ln2 - 1 + 



dx 



1 + 



1 

x 



(l + x) 



where v is the scale parameter in the MS-bar scheme. We have subtracted 
and added the logarithmic UV-divergence in the form J x (l + |x|) _1 . The 
momentum dependence of £ — S c is in the 0(T _1 ) terms. Although X — S c 
has been calculated in the Feynman gauge, its leading T contribution given 
by the first term in (|24] ) is in fact independent of the gauge parameter a. 
This can be seen by going back to (|17]), setting p = and subtracting the 
n = mode. We can easily add the contribution of the \(<j)*<j)) 2 interaction 
to the selfenergy. The complete one loop E — E c is obtained by g 2 — > g 2 + 4A 
in fl2§. 

Consider next the gauge field self energy tensor corresponding to the 
diagrams in Fig. 1, 



2y (2k + p)^2k + p) v 

ii + k 2 m 2 + k 2 )(m 2 + (k + p) 2 ) 



(25) 



The logarithmic divergence of each of the two contributing diagrams can- 
cels in the sum, which can be shown by using the Ward-Takahashi identity 
P/J^tiu{p) =0. In 1 + 1 dimensions this identity furthermore tells us there is 
only one independent component in H^, which we choose to be IIn. Sum- 
ming over the Matsubara frequencies and continuing p2 analytically to p° we 
get 



n 



ii 



g 



X 



2(1 + 2n) (2fci+pi) 



k [ 2uo 2ujuj' 
(l + n + n')(w + u/) (n - ri)(u - u') 



(to + Uj'f 



(u-u') 2 -pl 



(26) 



where now u 



+ k 2 



[u', n and n' are as in (pl|)). The dominant large 
T behavior of ELu is given by its classical approximation obtained with the 
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substitution 



9 2 T ^ 



(2fci+pi) 



(w + u/) 



A. 2 



(u; + u/) 2 - Po 



+ 



A2 



, (27) 



which is UV-finite and of order T. The next to leading behavior turns out 
to be of order T° (not InT): 



nn 



L1 n 
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, , + 0(T~ 

7T p( - P5 



(28) 



This contribution corresponds to loop momenta k of order T and it is there- 
fore the analogue of the hard thermal loop expression in 3+1 dimensions. 
Here in 1+1 dimensions it is subdominant to the classical Hq. All other dia- 
grams are superficially UV-convergent in two dimensions and for convergent 
diagrams the corrections to the approximation ([22]) are down by two powers 
T, as follows for example from the fact that they correspond to the n ^ 
terms in the summation over Matsubara frequencies. 

We conclude that in two dimensions the classical theory is a good ap- 
proximation to the quantum theory, for weak coupling. The classical theory 
is UV-finite. To minimize the difference between the classical and quan- 
tum theory we can match the classical mass parameter fi 2 of the Higgs field 
according to 



p 2 + S - S c 



-2 . 9' 
fi +- 



+ 4A 
2tt 



c — In 



T 



(29) 



where the bar on the right hand side indicates the MS-bar scheme. However, 
this matching is only of limited use since we cannot match the corresponding 
non-analytic subdominant terms in with a local classical action. 

An important point is now that we can take a limit of weak coupling and 
high temperature, such that the classical approximation becomes exact: 



A 



9 



£A, T 



oo. 



(30) 



The limit is such that X" and £ are kept fixed. Note that T' = AT/|/^| 3 , and 
AT and g 2 T are the combinations appearing in the dimensionally reduced 
theory. In this limit 
go to zero. In 



and (£8 



1 and the quantum corrections 
we have written \fi\ such that it also applies to the 



classical Higgs phase, for which fi < 0. Instead of /i we can of course use 
any convenient mass scale, e.g. the zero temperature Higgs mass. 
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If we do not take the above limit, the question arises, how accurate is the 
classical approximation. The above calculations are carried out for external 
frequencies and momenta which are small compared to the temperature. On 
the other hand, in real time (as opposed to frequency space) we also may 
expect a limited region of validity. The correction ( pgj) leads to propagation 
and situations might exist where it becomes substantial, depending on the 
physical quantity under study. In the application to the computation of the 
sphaleron rate, the Chern-Simons number is a zero momentum observable 
(explicitly so in Coulomb gauge), and the rate is should be dominated by 
low frequencies (cf. the cosine transform method in fT4||). 

The classical approximation is on a much better footing in 1+1 dimen- 
sions than in 3+1 dimensions, because it gives directly the dominant high 
temperature behavior of correlation functions, not only in the static case 
corresponding to dimensional reduction, but also in case of time dependence. 
(The classical partition function still suffers from the Rayleigh-Einstein- Jeans 
divergence.) In 4D, matching mass parameters between the original theory 
and the dimensionally reduced theory involves not only the classical combi- 
nations g 2 T and AT, but also g 2 T 2 and AT 2 . The latter combinations corre- 
spond to the well known hard thermal loops effects, which are the result of 
loop momenta of order of the temperature, and which reflect the quadratic 
divergences in the bare theory. Since g 2 T 2 dominates over g 2 T\ft\, the details 
of the matching formulas are important in four dimensions even in the weak 
coupling limit. In 2D, such matching can be ignored in the weak coupling 
limit (|g). 

4 Lattice regularization 

The classical partition function is still a functional integral over all <£>(x) and 
7r(x) and needs regularization. In the numerical simulations this is provided 
by a lattice, for which the action takes the form. 



S 




n=0 I 



+ \ (d - iA 0n )4> n \ 2 - \exp(-iaAi 



i a 



)<t>: 



Vt+l — 





(31) 
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where a is the lattice distance, L = Na and v 2 = —fi 2 /X. From now on we 
only consider the case /i 2 < 0. In the previous section we gave arguments that 
we can use this action in the classical approximation, provided we replace /z 2 
by an effective parameter /i 2 which is related to /2 2 in the MS-bar scheme by 
(|29f) . In the weak coupling-high temperature limit (|3~0D the difference between 
/i 2 and fi 2 may be neglected, which is what we shall do in the following. 
Accordingly, we drop the subscript c on /i 2 and v 2 . 

We use the Coulomb gauge to obtain a hamiltonian description. The 
gauge condition Ai n+1 — A\ n = states that A\ is independent of n and 
related to the Chern-Simons number C by, 



oA x = 2ttC/N . 

The Gauss constraint 

V A 4° = 4° 1 4- 4° 1 - 94° 

/ j L -^mn- rL n — m+1 ' - rl m— 1 ^■ ri - r 

n 

where j° = —j is the charge density, 

jo = -i<j)*D <f) + i(D (j)y< 

is solved explicitly as 



A u 



VIA, 



(32) 
(33) 

(34) 
(35) 

The canonical momenta 



under conditions with zero total charge Q = a J2 n Jr. 
are defined by n n = dL/d<f) n , etc. where S = J dt L. 

It will be convenient in the following to use \fi\ = y\v as the basic unit 
of mass and to use scaled quantities, indicated with a prime, 



(j) = V(j)' 



A, 



>VXA' 



Then 



S 
S' 



v 2 S', 

N-l 

dt' a' ^ 

n=0 

in 2 

2 



(36) 



(37) 



2£ 



(4 



+i 



/ |2 



exp(— ia'A' Xl x 



i / |2 / /2 



y n+l 



(38) 
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and the hamiltonian can be written as 

H = v 3 V\H' , (39) 

i+i ~ 0« I 



+«' (ki 2 - ±) + y E ^ A ^;° , (40) 

m,n=0 J 

/„ = ("^'nC +«)/«', (41) 



where Pgr, 7r' n and n' n are the canonical momenta conjugate to C, <p' n and 
4>'* n . The Poisson brackets are normalized to one, e.g. {(fi' m , n' n } = S mn . The 
classical partition function is given by 

Z = JdP c dC n(dvr'„ d< d0' n d^) exp (-^j (42) 

with T' defined by 

T = v 3 V\T' , (43) 



in accordance with (|30|) . The C integration is over the compact domain 
-N/2 < C < N/2 (mod N). We also define a scaled diffusion rate V by 

A(t) = ([C(t) - C(0)] 2 ) ->Tt= V't 1 , (44) 

for large times. 

It is easy to see that there are N classical ground states given by 

C = k = 0, 1, • • • , N - 1, <f/ = —= exp(i2nkn/N), 

\2 

and Pc = n' n = 0. Writing (p' n = p n exp(i9 n ) we can introduce a winding 
number of the scalar field by 

Y N-l 

w = — ^2 d6 n , d6 n = 6 n+ i - e n (mod 2vr) £ (-vr, 7r]. 

Z1T n=0 

Then also w = k in the ground states. Both C and w are defined modulo N. 
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For convenience we record various scaled quantities: 



E' 



rricf, = v V Am $ , vtla = vv Am a , 



m 



V2, 



m A 



(45) 
(46) 



which characterize the scaled theory given by the action S' or hamiltonian 
H' . We also have 



r 



r 



2 2 



(47) 



for the large temperature behavior (O) or ([14]). 



5 Numerical method 

It is straightforward to derive Hamilton's equations from the Coulomb gauge 
hamiltonian ( fl0|) . We use a second order Langevin procedure to generate 
initial conditions according to the canonical ensemble. In this procedure, 
also known as the Kramers equation method, Hamilton's equations are mod- 
ified by adding noise and friction terms to the equations for the canoni- 
cal momenta. If this is done straightforwardly local 'current conservation' 
= dojn + (Jn+i ~ 3n)/ a ~ dfij^ would be violated, since j° contains the 
momentum of the scalar field. Consequently the total charge Q would be 
nonzero, which is a problem for the evaluation of the Coulomb energy. To 
deal with these constraints (which are similar to the Gauss constraint in the 
temporal gauge) we follow ref. and use polar coordinates (f)' = pexp(i9). 
The random Langevin forces are then only applied to the gauge invariant 
variables p, p p , and not to the gauge variant 9, pg. See the appendix for 
a description of the algorithm. The random forces are also not applied to 
C and Pc either, which makes it possible to monitor thermalization of Pa- 
in practice the singularity of the polar coordinates can lead to problems at 
the origin p = 0. We dealt with this in the following way. The system is 
kicked at time intervals At' = h by the random forces. If the kick is large in 
some sense, the updated coordinates may get too close to the origin, causing 
large discretization errors. We avoid this by replacing the single step update 
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by a variable stepsize leapfrog algorithm for the integration of Hamilton's 
equations over the interval (£', t' + h). The Langevin stepsize is still given by 
h, however. 

The finite Langevin stepsize h introduces an error which we monitored 
by requiring the 'output temperatures' T' c and T', 



to be such that (3c,p differed less than 0.1 in absolute value from (3 = 1/T'. 
This criterion is based on the interpretation that small Langevin errors lead 
to an effective temperature, the output temperature. Then the expectation 
T' oc exp(— (3E' S ), with E' s = 2a/2/3 ~ 0.94 suggests that an absolute error in 
(3 of si 0.1, causes a relative error in V of ~ 10%. 

Having produced an independent configuration of p's and g's we took this 
as the initial condition for the real time integration of Hamilton's equations. 
For this we could use the original cartesian coordinates, as the condition of 
charge zero is easily satisfied by projecting regularly onto zero charge (this 
only involves changes of machine precision order). After real time integration 
the Langevin process was started again and the process was repeated untill 
sufficient statistics was obtained. 

Our configurations were actually microcanonical to some extent, because 
the Langevin processes were stopped (for historical reasons) at times that the 
total energy had its mean value. As expected, we found in a few checks that 
the true canonical ensemble gave the same results within errors. The friction 
parameter was taken / = 1, and we checked that results do not depend on 
/, as should be the case. 

An accurate fifth order Adams-Bashforth-Moulton predictor-corrector multi- 
step integration algorithm was used to keep the numerical drift in the total 
energy sufficiently small. Otherwise we found that the diffusion A turned 
out not to be linear in t' for large t' . 

6 Low temperature region 

The results that will be presented are for £ = g 2 /X = 4, for which m' A = 
2 = y/2m'^. Fig. H shows the behavior of (\4>'\ 2 ) as a function of (3. We 
see a minimum near (3 = 3 where the crossover takes place between the low 
temperature region ((3 > 3) and the high temperature region ((3 < 3). Recall 





(48) 
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Figure 2: 



1/12 



) versus (3. 



that even this 'low temperature region' can be interpreted as a region of high 
unsealed temperatures: T/m^ = v 2 /{\f2(3) — > oo as v 2 — > oo. 

For low temperatures the Chern-Simons diffusion A tends to be small 
and large times t were needed to get A > 1. For the evaluation of the rate 
we require A > 10, such that at least three sphaleron transitions take place 
(AC>3). 

In fig. |3] a comparison is made with the analytic semiclassical result (|12"D. 
The upper data are for a' = 0.32, the lower data for a' = 0.16 (at (3 = 10, 11 
for a' = 0.32 only). The system size is V = 16 (N = 50 and N = 100). We 
obtained the same results within errors for the smaller volume V = 10.28, 
in a check for a few (3 values. However, for sizes as low as V = 8 we did 
see clear deviations in the low temperature regime. The errors in the rate 
T are statistical and determined with the jackknife method from at least 
900 configurations. For (3 we used the output temperature ( fRf ), with errors 
obtained by the binning method. The input (3 is the nearest integer. 

We conclude from fig. ^| that the classical simulation is able to reproduce 
the semi-classical formula for /3^7. For the two lattice spacings a' = 0.32, 
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Figure 3: Results for lnF = ln(r/m|L) for £ = 4 and volume L' = 16. The 
diagonal solid line represents the analytic semiclassical result (12). The two 



solid lines extending from small (3 into the intermediate temperature region 
represent high temperature fits of the form F = cq + C2/3~ 2 0. The upper 
data are for a' = 0.32, the lower data for a' = 0.16. The figure also shows at 
very small (3 a fit of F = c + c 2 /3/3~ 2//3 to the a' = 16 data. 
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0.16 the classical rate in this temperature region is lattice spacing indepen- 
dent, within the errors. From the perturbative analysis in sect. 3 we see 
no reason to doubt that the limit of zero lattice spacing exists in the weak 
coupling region and that the resulting ratio r c j ass /r sem j_ c j ass approaches 
one as (3 — > oo. 

7 High temperature region 

Decreasing (3 we see in fig. |] that the rate is at first more or less temperature- 
independent, after which an apparent T 2 behavior sets in, which seems to 
confirm the expectation ([131) . However, there is substantial lattice spacing 
dependence in this temperature region and our earlier analysis lead to the 
conclusion that the coefficient c 2 in the fit F — co+c 2 T' 2 is proportional to a' 2 , 
such that the rate turned out temperature-independent in the limit of zero 
lattice spacing ||. Here we reconsider the situation at higher temperatures. 
We generated new data in the region 1.5 < T' < 5 (input (3 = 0.75, 0.65, 0.45, 
0.29, 0.23, 0.2), with volume L' = 16. At these temperatures this value of L' 
is amply sufficient for avoiding significant volume dependence. The lattice 
spacings where given by a' « 0.25, 0.23, 0.16, 0.11, 0.08 (a' = L'/N, N = 63, 
71, 100, 141, 200). At first we tried extrapolating to zero lattice spacing 
and seeing if the resulting rate showed a T 2 or T 2 / 3 behavior, but in this 
way the errors turned out too large to allow for a meaningful discrimination 
between these two alternatives. We therefore looked at the temperature 
dependence at fixed lattice spacing, with the following amusing outcome: the 
larger lattice spacings a' = 0.25, 0.23, favor T 2 behavior while the smaller 
spacings a' = 0.16, 0.11 and 0.08 favor T 2//3 behavior. Figs. || and [5] show 
two examples (as before T' in the plots is the output temperature). Fitting 
the form F = c + c 2 T' 2 (cf. (§7|)) led to x 2 /d.o.f. = 1.4, 0.42, 6.5, 4.8, 3.45, 
for a' = 0.25, 0.23, 0.16, 0.11, 0.08, respectively, while F = c + c 2/3 T /2 / 3 led 
to x 2 /d-o.f. = 28.9, 21.6, 2.5, 1.1, 1.2. Clearly, the T 2 / 3 form is favored for 
the three smaller lattice spacings. We extrapolated the resulting c 2 / 3 to zero 
lattice spacing, assuming a quadratic dependence on a' (cf. fig. [(J), with the 
result c 2 / 3 (a' = 0) = 0.00452(86), or 

k = 2c 2/3 = 0.0090(17), £ = 4. (49) 

It seems reasonable that the lattice spacing dependence falls like a 2 as a — » 0, 
since the equations of motion are based on an action with 0(a 2 ) discretization 
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5 10 15 20 1 1.5 2 2.5 

output T' a output T' 2/3 

Figure 4: Data showing T 2 behavior for relatively large lattice spacing, a' = 
0.23. Left: F versus T' 2 and a fit to C0 + C2T' 2 . Right: same data versus T l2 l z 
with the T' 2 fit of the left plot (dashed line) and a fit to c + c 2 /3T /2//3 (solid 
line). 

errors. This could be upset by divergences, such that odd powers of a appear, 
as in 3+1 dimensions, but these are not expected here since we have seen in 
sect. 3 that the classical theory is finite. 



8 Conclusion 

On the basis of perturbation theory we have argued that the effectively clas- 
sical approximation becomes exact in the limit (|30| ) of weak coupling and 



high temperature. Furthermore, for correlation functions the classical theory 
is finite, which means that we can in principle obtain continuum answers 
by extrapolating the numerical results to zero lattice spacing. This avoids 
having to determine an effective action with lattice spacing dependent coun- 
terterms, as we previously deemed necessary || . In the limit ([It]) the purely 
classical action is sufficient to leading order. 

In the resulting classical model we can still distinguish a low and high 
temperature regions separated by a crossover domain, as illustrated by the 
behavior of the rate in Fig. [| In the low temperature region the rate is well 
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5 10 15 20 1.5 2 2.5 

output T' E output T' 2/3 

Figure 5: Data illustrating the T 2 / 3 behavior at small lattice spacings, a' = 
0.16. Right: F versus T /2//3 and a fit to Cq + C2/?,T' 2 ^. Left: same data versus 
T' 2 with the T' 2 / 3 fit of the right plot (dashed line), and a fit to c + c 2 T' 2 
(solid line). 



described by the analytic form based on the sphaleron fllQ |, as observed earlier 
@, |7|, H . The lattice spacing dependence is small in this region, of the order of 
the error bars (~ 10%). Upon crossing over to the high temperature region 
the lattice spacing dependence grows, and leads to a tricky phenomenon: 
apparently Y/L oc T 2 at intermediate lattice spacings, but this turns into 
a T 2 / 3 behavior in the continuum limit. We obtained the estimate Y/L = 
0.0090(17)(AT) 2 / 3 , T -> oo, for £ = g 2 /X = 4. 

As mentioned in the Introduction, arguments for a T 2//3 behavior were 
previously given in ref. f|. It was shown there that the classical rate could 
be written in terms of a function g(x,y), as Y/L = T 2 ^ 3 g(a 3 T,v 3 /T), using 
units in which A = 1, suppressing £ dependence. For high temperatures T 2 / 3 
behavior then followed in the continuum limit provided g(0, 0) 7^ 0, 00 ex- 
isted. For the case v 2 = 0, numerical data supported a nontrivial g(0,0). 
However, for v 2 = 1 (which corresponds to our scaled units fl56p) the results 
for the rate were not compared to T 2 / 3 behavior with much precision. The 
data looked compatible with our's at the time |J, but since it was taken at 
different £ (£ = 10) a detailed comparison was not meaningful. Our analytic 
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results in Sect. 3 support the existence of a continuum limit of g(a 3 T, v 3 /T), 
but notice how the combination a 3 T implies non-commutativity of the limits 
a — > and T — > oo. Such non-uniform behavior complicated our numerical 
analysis indeed. 
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A Simulation algorithm using Kramers equa- 
tion 

This appendix we describe the second order Langevin algorithm, closely re- 
lated to 'Kramers equation'. Consider the class of systems described by the 
lagrangian 

L=\g af i(q)q a q f) -V(q), (50) 

where g a p plays the role of a metric in coordinate space. We shall assume 
that g a p can be expressed in the form 

g a p = e ak ef3 k , a,P,k — l,...M, (51) 

where M is the number of g's or p's. The abelian Higgs model in polar 
coordinates can be described this way. It is useful to have in mind a simple 
two dimensional illustration given by 

L = 4>*4>-V{<t>) (52) 
= f + fp-vfag), <p = p e\ (53) 

for which 

en = y/2, e 22 = V2 p, e 12 = e 2 i = 0. (54) 
The hamiltonian is given by 

H = \g aP p a p P + V, p a = g afS qp, (55) 
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where g al3 is the inverse of g a p, g a ^g^ = 5". 

Discretizing time with stepsize h, t = nh {n = 0, 1, . . .), we can derive a 
leapfrog algorithm by requiring the action 

S = 5>«m(Ci - q a n ) - hH(p n , q n )\ (56) 

n 

to be stationary under variations of the canonical variables, 

dH(p n+1 ,q n+1 ) 

Pan+1 = Pan ~ fl ^— . (58) 

We modify these equations into second order Langevin form by adding fric- 
tion and noise terms to the momentum equations, 

Pan+1 = Pan ~ tl — hF£(q n+1 )p 3n 

VQn+l 



+V2hN ak (q n+ i)u kn , (59) 
where the v kn are normalized gaussian random numbers, 

\yhn) = 0, (Vkn^k'n 1 ) = $kk'$nn'- (60) 

By a suitable choice of the random forces and friction the system may equili- 
brate (at large times and in the limit of zero stepsize) according to the Boltz- 
mann distribution exp(—H/T). For this it is necessary that the Boltzmann 
distribution is a fixed point of the associated Fokker-Planck-type equation, 
which is called Kramers equation [31]. This 'stability condition' leads to 
relations between F and N, which will now be derived. 

By the usual arguments one can derive the evolution equation for the 
probability distribution W t (p,q) at time t ('Kramers equation') 

w t+h - w t (_^EJL + ^EJL (en 

h \ dp a dq a dq a dp a 

^ + f ^ + M ^Sf, + o(ft) ) w " (62) 



where 



M a/3 = N ak Nfs k . (63) 
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The left hand side approaches W in the limit of zero stepsize and the first 
two terms on the right hand side represent the usual Liouville flow, for which 
exp(— H/T) is a fixed point. For the other terms exp(— H/T) is also a fixed 
point provided that 

»-< :-x»*k+"* {M^-^Ba)- (64) 

Inserting the explicit form of H and comparing powers of the momenta leads 
to the relations between F and N, 

K = \M aP g^, (65) 
I(F 7 V^ + F 7 V") = ^M^g*. (66) 



This is satisfied by 



FI = \m^, (67) 



which gives F in terms of M and the metric, 

Ergodicity and stability should guarantee that the evolution approaches 
the canonical distribution for t — > oo. This can be understood intuitively by 
writing 

W t = e~ H/2T W t . (68) 
The distribution W t evolves according to 

w M -w t + (d H a m e\ Wt = _ m , m 



with 



h \dp a dq a dq a dp a/ 



+ )M aP — + . (70) 



V d Pa 2T J ^\d Pf} 2T 

The Liouville operator on the left hand side conserves probability. On the 
other hand, since M is evidently a positive matrix, the differential operator 
P is positive definite. Expanding W t in terms of the eigenmodes of P, the 
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modes with nonzero eigenvalues are expected to die out during the time 
evolution, such that W t approaches the zero mode exp(— if/2T), 

13 +<^W = 0, (71) 



\dp„ 2T 

which is of course also a zero mode of the Liouville operator. The original 
Wt then approaches the Boltzmann distribution with temperature T. 
A natural choice for N is given by 



N ak = ^fTe ak , M a/3 = fT 9aP , FZ = f6*, (72) 

with some arbitrary friction coefficient / > 0. However, in our application to 
the abelian Higgs model we want to add noise only to the radial momenta p p 
and not to the angular momenta pg. Therefore, N and M contain a projector 
onto the radial variables, 



N a k = y fT e a iiTi k , Tr k iir lm = 7r fcm , (73) 

where 7r is the projector. In the above two dimensional model this is illus- 
trated by 

7Tn = 1, 7Ti2 = 7T 2 l = 7T 2 2 = 0, (74) 

and consequently 



N n = ^2fT, Mn = 2/T, F* = f, (75) 

with the other components vanishing. These N and F are independent of p 
and 9. 
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